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Abstract. In the last decade ground-based Imaging Atmospheric Cherenkov Telescopes have 
discovered roughly 30 pulsar wind nebulae at energies above 100 GeV. We present first results 
from a leptonic emission code that models the spectral energy density of a pulsar wind nebula by 
solving the Fokker-Planc.k transport equation and calculating inverse Compton and synchrotron 
emissivities. Although models such as these have been developed before, most of them model the 
geometry of a pulsar wind nebula as that of a single sphere. We have created a time-dependent, 
multi-zone model to investigate changes in the particle spectrum as the particles diffuse through 
the pulsar wind nebula, as well as predict the radiation spectrum at different positions in the 
nebula. We calibrate our new model against a more basic previous model and fit the observed 
spectrum of GO.9+0.1, incorporating data from the High Energy Stereoscopic System as well 
as radio and X-ray experiments. 


1. Introduction 

In the last decade ground-based Imaging Atmospheric Cherenkov Telescopes (IACT S ) have 
discovered almost 150 very-high-energy (VHE; E > 100 GeV) 7 -ray sources. Roughly 30 of 
these are confirmed pulsar wind nebulae (PWNe), while other source classes include supernova 
remnants, active galactic nuclei, or unidentified sources Q A subset of the unidentified sources 
may eventually turn out to be PWNe. H.E.S.S. is currently the world’s largest IACT and is 
located near the Gamsberg mountain in Namibia. H.E.S.S. consists of five telescopes which were 
built in two phases. In December 2003 Phase I consisting of four 12-m telescopes was completed. 
The fifth 28-rn H.E.S.S. II telescope has been operational since July 2012, making H.E.S.S. the 
largest and most sensitive ground-based 7 -ray telescope in the world. 

A PWN is a bubble of shocked relativistic particles, produced when a pulsar’s relativistic 
wind interacts with its environment [Tj. The conversion of a fast-moving particle wind into 
electromagnetic radiation happens by two main processes. The leptons in the particle wind 
interact with the magnetic field of the nebula and this causes synchrotron radiation (SR) up 
to several keV. The second process is when the high-energy leptons interact with low-energy 
photons through inverse Compton scattering (IC), boosting the photon energies up to E > 100 
GeV. Due to these two effects the radio, X-ray, and VHE 7 -ray emissions are tightly linked, as 
all three emerge from the same lepton population. The following are typical characteristics of a 
PWN 0: 
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(i) A filled-centre morphology which is brightest at the centre and dimming in all directions 
towards the edges; 

(ii) A flat radio photon spectrum with an index between 0 and -0.3; 

(iii) A well organised internal magnetic field structure; 

(iv) A high level of linear polarisation at high radio frequencies; 

(v) Evidence of particle re-acceleration; 

(vi) Evidence of synchrotron cooling which means that the size of the PWN decreases with 
increasing energy. 

There are many unanswered questions in PWN physics. Surprisingly, it was noted [3j that 
the measured 7 -ray luminosity (1-10 TeV) of the PWNe does not correlate with the spin-down 
energy of their embedded pulsars. On the other hand, it was also found [1] that the A'-ray 
luminosity is correlated with the pulsar spin-down energy. Furthermore, it is currently unknown 
whether there is any correlation between the TeV surface brightness of the PWNe and the spin- 
down energy of their pulsars. Due to this reason it is necessary to create a spatially dependent 
model to calculate the spectral energy density (SED) from the PWN. The spatial dependence 
will yield the flux as a function of the radius. We calibrate our model against a more basic 
previous model and fit the spectrum of a known PWN. We describe our model in Section [2j 
and its calibration and spectral fitting to GO.9+0.1 in Section [3} We discuss our conclusions in 
Section [4] 


2. PWN Model 

2.1. Transport equation and injection spectrum 

We model the transport of charged particles in a PWN by solving a Fokker-Planck-type equation 
similar to the Parker equation J3] . This equation includes diffusion, convection, energy losses, 
as well as a particle source. If we neglect spatial convection, we are left with (5] 
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with Q the particle injection spectrum, N e the particle spectrum, E e the particle energy loss 
rate, and k the diffusion tensor. We used a broken power law for the particle injection spectrum 
following [ 2 ], 


Q(E e ,t) 
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Here Qq is the normalisation constant, Eq the break energy, E e the lepton energy and a± and 
a 2 the spectral indices. To obtain Q 0 we use a spin-down luminosity L(t) = Lq/ (1 + t/r, o ) 2 of 
the pulsar [ 8 ], with to the initial spin-down time scale of the pulsar, and Lq the initial spin-down 
luminosity. We normalise Q as follows: 
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with e the conversion efficiency of the spin-down luminosity to particle power. 




2.2. Radiative energy losses in the PWN 

Particle energy is dissipated from the system due to radiation. We incorporated SR and IC 
scattering, similar to calculations done by [B] in their globular cluster model. SR losses are 
given by [9] 

( dE e\ U T C E 2 B 2 

V dt ) SR 6 vr (m e c 2 ) 2 ’ 

with a t the Thompson cross section and B the PWN magnetic field. The IC scattering energy 
loss rate is given by 
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with gic = 2ir e 4 c, e the soft photon energy, n £tP the blackbody photon number density, 
T p the photon temperature of component p, E 1 the TeV up-scattered photon energy, and 
((E e , Ej,e) = £oC(Re, R 7 , e) the collision rate with £0 = 2-Ke i Eoc/ eE 2 , and ( given by 
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Here, £ 0 = "i e c 2 , f(q,g 0 ) = 2 glng + (1 - ?)(1 + (2 + go)?), ? = RqR 7 /(4e£ , e (R e - R 7 )), and 
< 7 o(£, R 7 ) = 2eE 1 /Eq. The two radiation loss rates can be added to find E e used in 0. 


2.3. Calculation of the particle spectrum 

If we use spherical coordinates and assume spherical symmetry, i.e. = 0 and = (1, the only 
change in particle (lepton) spectrum will be in the radial direction, so V 2 = \ (J^ ( r 2 yp)) • 
For the diffusion scalar coefficient k we consider Bohm diffusion 
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with kb = c/3e, and c and e denote the speed of light and the elementary charge. We can now 
discretise Q: 
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with z the time index, j the energy index, k the radial index, /3 = 2ndt/Ar 2 , 7 = 2ndt/(rAr) 
with Ar the size in the spatial bins. Also r a = AEi+i/AEi with 


Z -(aE i+1 -AeJ (r, (9) 

We initially approached the discretisation process by using a simple Euler method. It soon 
became clear that this method was not stable. We then decided to use a DuFort-Frankel scheme 
to solve 0. In solving this equation, we calculate the electron spectrum of the PWN due to 










the injected particles from the embedded pulsar, taking into account their diffusion through the 
PWN and the IC and SR energy losses. We use the following parametrised form for B(t) [2] 


B{t) 


1 + 



( 10 ) 


with B 0 the birth magnetic field, t the age of the PWN, and as the magnetic field parameter. 
We can calculate Lq and Bq by using the present-day luminosity and the present-day magnetic 

field together with the age of the PWN. We limit the particle energy using E max = § yj c{i+^ ) [2j> 
with a the ratio of electromagnetic to particle luminosity. Particles with E e > B max are assumed 
to escape. 

To solve ([ 8 ]) numerically, our multi-zone model divides the PWN into shells. The particles 
are injected into zone one close to the centre and are allowed to diffuse through the different 
zones. We assumed as initial condition, that all zones were devoid of any particles. We used a 
reflective boundary condition at the inner boundary of the PWN, and on the outer rim we set 
N e equal to zero to allow the electrons to escape from the PWN. 


2-4- Calculation of radiation spectrum 

The time-dependent photon spectrum of each zone can now be calculated, given the electron 
spectrum solved for each zone. For IC we have [ 6 ] 
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where A = 47 rd 2 , d the distance to the source, and A f e is the number of electrons per energy in a 
spherical shell around r. We consider three components of target photons, cosmic background 
radiation (CMB) with T\ = 2.76 K and u\ = 0.23 eV/cm 3 , Galactic background infra-red 
photons with T 2 = 35 K and U 2 = 0.5 eV/cm 3 , and starlight with X 3 = 4 500 K and = 50 
eV/cm 3 . For SR we have 
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Table 1. Values of model parameters. 


Model Parameter 

Symbol 

Value 

Braking index 

n 

3 

.B-field parameter 

as 

0.5 

Present day R-field 

B(T) 

40.0 p,G 

Conversion efficiency 

e 

0.6 

Age 

T 

1900 yr 

Initial period of the pulsar 

Po 

43 ms 

Current period of the pulsar 

P 

52 ms 

Pulsar spin-down rate 

P 

1.5557 x 10" 13 

Distance 

d 

8.5 kpc 

Qindex 1 

OL 1 

-1.0 

Qindex 2 

<7 2 

-2.6 

Maximum energy parameter 

a 

0.2 

Soft photon component 1 

T\ and u\ 

2.76 K, 0.23 eV/cm 3 

Soft photon component 2 

T 2 and u 2 

35 K, 0.5 eV/cm 3 

Soft photon component 3 

T 3 and U 3 

4500 K, 50 eV/cm 3 
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Figure 1. Calibration of our model against the model of [2j. The radio data are from mi, the 
X-ray data are from m , and the 7 -ray data are from |13j . 


with u cr the critical frequency (with pitch angle 9) 

Vcr(Ee,0,r) = ^E*B(r) sin. 0 (13) 

and k given by 

poo 

k(x) = x I K 5/3 (y)dy, (14) 

J X 

where K 5 / 3 a modified Bessel function of the second kind of order 5/3. 

2.5. Calculation of the line of sight flux 

In order to do the line of sight (LOS) integration of the flux from the PWN we need to use 
the flux per unit volume calculated in the previous section and multiply it with the volume in 
a particular LOS as viewed from Earth. Since our model is spherically symmetric and since 
the source is very far from Earth, we used cylinders intersecting the spherical zones. We chose 
the cylinders and the spheres to have the same radii. We then used the intersection volume 
to calculate the flux in a particular LOS. We will present the radial dependence of the flux 
elsewhere. 

3. Model calibration and spectral modelling of GO.9+0.1 

To calibrate our model we compared our results to those of [2j as shown in figure [I] Pulsar 
J1747—2809 has now been discovered in PWN GO.9+0.1 [TO] with P = 52 ms and P = 
1.5557 x 10' 13 . We used these values to calculate tq = Pq/2Pq assuming 1+ = 43 ms and 
no decay of the pulsar B-field. In the figure, the radio data are from DU, the X-ray data are 
from [ 12 ] . and the 7 -ray data are from [13]. We list our model parameters in table [I] 










4. Conclusions 

We created a time-dependent, multi-zone SED model to model the radiation spectrum observed 
from PWNe. We calibrated this model using a previous model. We then fitted the data of 
GO.9+0.1. The new model does not exactly reproduce [2], but the results are quite close. The 
new SR component is lower due to the fact that the old model did not consider IC losses while 
performing the particle-transport. We will next perform a population study to probe a potential 
relationship between the TeV surface brightness and the spin-down luminosity of the embedded 
pulsar. We also intend to investigate the spatial dependence of the different fluxes. 
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